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We numerically investigate the transport of a suspended overdamped Brownian particle which is 
driven through a two-dimensional rectangular array of circular obstacles with finite radius. Two 
limiting cases are considered in detail, namely, when the constant drive is parallel to the principal or 
the diagonal array axes. This corresponds to studying the Brownian transport in periodic channels 
with reflecting walls of different topologies. The mobility and diffusivity of the transported parti- 
cles in such channels are determined as functions of the drive and the array geometric parameters. 
Prominent transport features, like negative differential mobilities, excess diffusion peaks, and un- 
conventional asymptotic behaviors, are explained in terms of two distinct lengths, the size of single 
obstacles (trapping length) and the lattice constant of the array (local correlation length). Local 
correlation effects are further analyzed by continuously rotating the drive between the two limiting 
orientations. 



PACS numbers: 05.40.-a,05.60.Cd,51.20.-|-d 

I. INTRODUCTION 

The effective control of mass and charge transport in 
artificial micro- and nanostructures requires a deep un- 
derstanding of the diffusive mechanisms involving small 
objects in confined geometries. Such situations are typi- 
cally encountered when studying the transport of parti- 
cles in, e.g.: biological cells pj and zeolites [2], catalytic 
reactions occurring either on templates or in porous me- 
dia [J, chromatography or, more generally, separation 
techniques of size-dispersed particles on micro- or even 
nanoscales [4]. In many respects these transport phe- 
nomena can be regarded as diverse manifestations of ge- 
ometrically constrained Brownian dynamics in one (ID) 
or higher dimensions [S]- 

In this paper we focus on the dc driven transport of 
pointlike Brownian particles in two dimensional (2D) ar- 
rays. Restricting the volume of the phase space available 
to the diffusing particles by means of confining bound- 
aries, or obstacles, causes remarkable entropic effects 6 - 
117] . In quasi-lD geometries (narrow channels), driven 
transport of charged particles across bottlenecks, such 
as ion transport through artificial nanopores or biolog- 
ical channels, represents an ubiquitous situation, where 
diffusion is effectively controlled by entropic barriers [S]. 
Similarly, the operation of artificial Brownian motors and 
molecular machines [9", '10] also results from the interplay 
of thermal diffusion and the pinning action by both en- 
ergetic and entropic barriers. 

Higher dimensional geometries have recently attracted 
broad interest in the context of separation of macro- 
molecules [TS], like colloids [12], DNA fragments [2011^ . 
or even magnetic vortices p21 - B5] , because of the occur- 
rence of induced transverse drifts, which separate differ- 
ent objects depending on their bulk diffusivity [26M8] . 
Such systems are essentially 2D arrays of impenetrable 



obstacles traversed by diffusing particles subjected to ex- 
ternal gradients. An accurate modeling of real experi- 
ments, like those cited here and many more, would re- 
quire incorporating nontrivial effects due to particle in- 
teractions ^22, ,27j, fluidics [Ml US], chaos [5^, inertia 
[S5] , excluded volume HI] , particle shape |H2] , spa- 
tial asymmetries and disorder [9l [10] , to mention but a 
few. In other contexts the underlying planar constrained 
geometries have been modeled also by arrays of traps 
[551 151) , or egg-carton potentials [33 [3S] . 

In order to attempt a first quantitative characteriza- 
tion of forced transport across a 2D array, we used a 
simplified model where a single Brownian particle of neg- 
ligible size is suspended in an unmovable interstitial fluid 
at fixed temperature (dilute suspension). The particle is 
free to diffuse in the connected space delimited by cir- 
cular reflecting obstacles of finite radius, arranged in a 
rectangular lattice. The particle is subjected to ther- 
mal fluctuations and large viscous damping (as is often 
the case for biomolecules, colloids and magnetic vortices 
[21 [TU]), and a homogeneous constant force (Fig. [ij. Such 
a dc drive is applied from the outside by coupling the 
particle to an external field (for instance, by assuming 
that the particle carries a dielectric or magnetic dipole, 
or a magnetic flux) , without inducing drag effects on the 
suspension fluid. To simplify our model even further, 
we neglected any local spatial asymmetry, such as one 
obtains by lining up non-circular obstacles in symmet- 
ric lattices or by arranging circular obstacles in suitable 
anisotropic arrays [3 HU] . Under such conditions the dy- 
namics of the driven particle is fully described by an 
overdamped 2D Langevin equation, which we encoded 
in our numerical simulation algorithm, or by the corre- 
sponding 2D Fokker-Planck equation, which lends itself 
to a more systematic analytical treatment [5]. All other 
effects, including particle-particle and particle-obstacle 
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FIG. 1: (Color online) Sketch of a square array, ~ Ly, of circular obstacles with radius tq. (a) The Brownian particle is 
driven by a dc force F oriented at an angle </!> with the horizontal axis. Transport through the array can be reduced to transport 
along corrugated channels of type I for F^, = 0, and type II for F^^, (f> = 7r/4, with appropriate compartment sizes. Channels 
I and II are sketched in (b) and (c), respectively; y^" and x^" denote the width and the periodicity of the channel. The vertical 
distance between two opposite obstacles is 2A^''' = y^" — 2ro. The trajectory drawn in (c) is the limiting noiseless trajectory 
for F^^ — >■ oo pointing to the right. The central horizontal lines in (b) and (c) are the channel symmetry axes, which suggest 
to further reduce the analysis of the problem to half the sketched channels (see Sec. [ll| . 



interactions, hydrodynamic corrections, chaotic and in- 
ertial dynamical terms, particle size and shape, are thus 
ignored. 

The properties of Brownian transport in the stationary 
regime are well quantified by the particle mobility and 
diffusivity as functions of the external drive (magnitude 
and orientation) and of the array parameters (obstacle 
radius and lattice constants). Thanks to the simplifica- 
tions assumed in the present model, we detected promi- 
nent transport features, like negative differential mobil- 
ities, excess diffusion peaks, and unconventional asymp- 
totic behaviors, which went unnoticed in earlier reports, 
experimental and computational, alike. We explain the 
properties of Brownian transport in 2D arrays in terms 
of two distinct depinning mechanisms: (i) trapping by a 
single obstacle, which the particle can overcome only by 
diffusing a transverse distance of the order of the obsta- 
cle size, while being driven the same longitudinal length; 
(ii) correlated collisions against obstacle rows at an angle 
with the external force, which tend to collimate the parti- 
cle trajectories to form stream lines connecting and flow- 
ing around the obstacles. Such mechanisms correspond 
to two different length scales, respectively: the obstacle 
radius (independent of the array geometry) and the ef- 
fective obstacle spacing (resulting from the combination 
of array geometry and drive orientation). 

This paper is organized as follows. In Sec.|II]we intro- 
duce the Langevin equation formalism employed in our 
simulation code. We show how for special orientations 
of the drive the problem can be reduced to the trans- 
port in corrugated ID channels. Mobility and diffusivity 



obstacle radius, and the array lattice constants for two 
such corrugated channels corresponding, respectively, to 
orienting the drive along the principal (Sec. Ill A) and 



the diagonal axes of a square array (Sec.|IIIB). Mobility 



and diffusivity dependence on the system parameters are 
interpreted in Sec. |IV| in terms of the two aforementioned 
depinning mechanisms. In Sec.|V] local correlation effects 
due to the array geometry are analyzed in more detail by 
continuously rotating the drive between the two limiting 
orientations of Sec. |III[ Finally, we summarize our results 
in Sec. IVD 



II. MODEL 

Let us consider an overdamped Brownian particle of 
unit mass diffusing in a suspension fluid contained in a 2D 
rectangular array, L^x Ly, of reflecting circular obstacles 
of radius rg, as illustrated in Fig. ^a). The particle is 
subjected to a homogeneous force F oriented at an angle 
(f> with the horizontal axis x. The overdamped dynamics 
of the particle is modeled by the 2D Langevin equation. 



df 
di 



= F 



(1) 



data are plotted in Sec. Ill as functions of the drive, the 



where f — {x,y) and £^{t) — {^x{t),^y{t)) are zero- 
mean, white Gaussian noises with autocorrelation func- 
tions {ii{t)^j{f)) = 2S,jS{t - t'), with i,j = x,y. The 
noise strength Dq corresponds to the particle free or bulk 
diffusivity in the absence of geometric restrictions and is 
proportional to the temperature of the suspension fluid. 
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FIG. 2: (Color online) Driven transport in channel I with 
XL = vl ~ ^, Dq = 0.03, and different ro: particle mobility 
fi (a) and scaled diffusivity D/Dq (b) vs. the scaled force 
F/Dq. Inset: mean exit time through four pores, r4, and two 
opposite pores, T2. The numerical data are compared with 
the analytical predictions in Eq. ([8|. 



We numerically integrated Eq. ([T]) by a Milstein algo- 
rithm [37] . The stochastic averages reported in the forth- 
coming sections were obtained as ensemble averages over 
10^ trajectories with random initial conditions; transient 
effects were estimated and subtracted. To check the re- 
liability of our numerical simulations we solved the cor- 
responding Fokker-Planck equation by means of a finite- 
element algorithm and obtained a good comparison with 
the Langevin method. 

The driven particle current across the array strongly 
depends on the force orientation. The problem cannot be 



reduced to a ID problem, unless F is oriented along a lat- 
tice axis, i.e., (j) = 4)n,m with ipn.m = aictan(nLy/mLx) 
(n,m S Z and m ^ 0). Due to its spatial symmetry, 
the array can then be regarded as consisting of iden- 
tical periodic channels parallel to F and separated by 
reflecting walls. This allows to interpret the stationary 
transport properties of the array in terms of transport 
mechanisms through compartmentalized narrow (quasi 
ID) channels. In Fig. [l] we illustrate the two limit- 
ing cases with = and (j) = ■'r/4 for a square lattice, 
Lx = Ly. For a rectangular lattice, i.e. Lx ^ Ly, and 
(/> = 0, the reduced channels are horizontal and have rect- 
angular compartments, a;^ x y\ with = Lx and = 
Ly. For (j) = 01 1 ~ aictan Ly/Lx, the reduced chan- 
nels are diagonal and their compartments have width 

fl = 2LxLyl^Ll + Ll and period = ^L^+Ll. 

For a square lattice, the latter channels coincide with the 
channel II of Fig. life), where y" = x'^l = LxV^- Note 
that for symmetry reasons, the discussion of channels I 
and II can be further simplified by halving the channel 
width, see Figs. [T|b) and (c). 

Motivated by these observations we investigated Brow- 
nian transport in two categories of reduced periodic chan- 
nels, conventionally directed along the x axis (and par- 
allel to F). They correspond to cutting respectively the 
channel I of Fig. [ijb) and the channel II of Fig. [IJc) 
along the central dotted line and then varying the com- 
partment length, XL, at will (note that here and in what 
follows the superscripts are omitted) . Channels I are thus 
characterized by a corrugated and a smooth wall, a ge- 
ometry that lets straight particle trajectories through, 
no matter what the radius of the circular obstructions. 
Channels II have equally corrugated walls, though shifted 
by half a period; for ro > yL/4 the driven particle can 
cross the channel only by circumventing the obstructions 
on either walls. As a consequence, one expects distinct 
transport properties for these two channel geometries 
[55ti43| . Note that tuning the compartment parameters 
corresponds to investigating rectangular (channels I) and 
face-centered rectangular arrays in 2D (channels II) with 
F parallel to the principal axes. 

Two transport quantifiers that best illustrate the dif- 
ferent properties of channels I and II are the nonlinear 
mobility and the effective diffusion coefficient. We char- 
acterize the response of a Brownian particle dc-driven 
along the channel axis by computing its mobility /i, 

M(F) = {xiF))/F^ (2) 

where {x{F)) = \imt^ac[{x{t)) — x{0)]/t, and its diffusiv- 
ity D, 

D{F)/Do = lim [(x2(i)) - {x{t))^]/2Dot, (3) 

r— >-oo 

both as functions of F for different channel geometries. 
In the absence of external drives, Einstein's relation |38) . 

^,o = m = D{0)/Do, (4) 



4 




FIG. 3: (Color online) Driven transport in a channel I: par- 
ticle mobility fi (a) and scaled diffusivity D/Do (b) vs. the 
scaled force F/Do for different xl- Other simulation param- 
eters are: j/l — 1, Do — 1, and ro — 0.4. 



establishes the dependence of the transport quantifiers on 
the channel geometry and temperature under equilibrium 
conditions. 

For a more compact presentation of our numerical 
data, we remark that Eq. ( T]) can be conveniently rewrit- 
ten in terms of the rescaled units t -> tDo and F — > 
F/Dq. a straightforward dimensional argument shows 
that both the particle mobility, Eq. ([2]), and its diffusiv- 
ity in units of Dq, Eq. ([3|, are functions of F/Dq, only, 
for any given channel geometry. 



III. CHANNEL TRANSPORT 

Contrary to smoothly corrugated channels, also called 
entropic channels [S], introduced first in Ref. [5] and 
further investigated in Refs. [TTHT5] . strongly compart- 
mentalized channels with narrow sharp bottlenecks, or 
pores, cannot be analyzed in terms of an effective ID ki- 
netic process directed along their axis (41j . Accordingly, 
driven transport in such strongly constrained geometries 
exhibits distinct features, which cannot be reduced to 
known properties of Brownian motion in ID periodic sys- 
tems [H1I3H]- 

We outline here the main differences between transport 
in entropic and sharply compartmentalized channels. In 
entropic channels fJ.{F) increases from a relatively small 
value for F = 0, of Eq. (|4L up to the free-particle 
limit, — 1, for — > oo (l41 115| . On the contrary, 
in coaxial compartmentalized channels n{F) decreases 
monotonically with increasing F towards a geometry- 
dependent asymptotic value, //oo, equal to the ratio of the 
pore to the channel cross-section ^IT . Such an asymptote 
vanishes for eccentric compartmentalized channels with 
off-axis, non-overlapping pores [15] . 

Significant differences have also been reported for the 
diffusivity. For entropic channels with smooth pores, 
the function D{F) approaches the free-diffusion limit for 
F — T' oo, D{oo) = Do, after going through an excess dif- 
fusion peak centered at an intermediate (temperature de- 
pendent [MKI^) value of the drive. Such a peak signals 
the depinning of the particle from the entropic barrier 
array [44,. In coaxial compartmentalized channels, in- 
stead, D{F) exhibits a distinct quadratic dependence on 
F [40l|42], reminiscent of Taylor's diffusion in hydrody- 
namics. This observation suggests that the particle never 
frees itself from the geometric constriction of the com- 
partment pores, no matter how strong F. The quadratic 
divergence D{F) can be cut off for suitably large F by 
alternately shifting the pores off-axis (eccentric channels 

m)- 

Driven transport in channels I and II of Fig. [l] is 
characterized by the superposition of, or sometimes the 
competition between, the properties observed for the two 
opposite compartmentalization geometries considered in 
the earlier literature. We remind here that, due to the 
mirror symmetry of the obstacle geometries of Figs. [l|b) 
and (c), our analysis refers to the irreducible channels I 
and II obtained by cutting the channels sketched there 
along their axis (dotted line) and then varying their spa- 
tial periodicity. Therefore, the irreducible channels dis- 
cussed below have effective width 2/_l/2 and transverse 
bottleneck width A = — tq. 



A. Channels I 

Channels I corresponding to square arrays, xl — Ul, 
exhibit transport properties apparently not much differ- 
ent from the entropic channels. In particular, the mobil- 
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FIG. 4; (Color online) Driven transport in a channel II with 
XL ~ 2, ul ~ 1, Do = 0.03, and different ro: particle mobility 
fi (a) and scaled diffusivity D/Dq (b) vs. the scaled force 
F/Dq. Dotted lines represent ^oc, Eq. Q, in (a) and Doo, 



Eq. (10 1, in (b) 



ity is an increasing function of F; the ^i{F) curves jump 
from fiQ to /Zoo = 1 sharper and sharper as the radius tq 
of the obstacles is increased [Fig. [2ja)]. The diffusivity 
curves develop an excess diffusion peak in correspondence 
with the mobility surge; the height of these peak weakly 
depend on tq. 

New features start to emerge when we consider chan- 
nels with rectangular compartments, i.e. XL > Ul- These 
become apparent in Fig. [3] where we display data for 
simulated transport in channels I for different period 



lengths XL. The occurrence of two distinct diffusion 
mechanisms is proved by the two-step increase of the 
mobility and the corresponding double diffusivity peaks. 
Such peaks merge a,s xl tends to ul- 



B. Channels II 

In Fig. [4] we plot the curves ij{F) and D{F) for 
channels II with comparable compartment dimensions, 
XL and Ul, and different obstacle radius, rg. The mo- 
bility function, ^J-{F), develops a non-monotonic depen- 
dence on F, which becomes prominent for narrow bottle- 
neck widths, A — yL/2 — tq, but vanishes altogether for 
A > Vl/'^ [Fig. |4f^a)]. Negative differential mobility is a 
peculiar feature of all sharp compartmentalized channels 
[H] and eccentric channels, in particular [43 . Indeed, 
each compartment of the reduced channels II of Fig. [2][c) 
has two off-axis bottlenecks, located alternately against 
its top and bottom wall. The subsequent mobility surge 
toward a horizontal asymptote with < 1 indicates 
that on increasing F the particle grows more sensitive to 
the pinning action of the obstacles, which it overcomes 
initially by mere Brownian diffusion and eventually with 
the assistance of the drive itself. 

Correspondingly, such change in the depinning mech- 
anism is signaled by a conspicuous diffusion peak [Fig. 
|4]Jb)]. In comparison to Fig. [2jb), here the D{F) peaks 
are more than one order of magnitude larger than Dq 
and their position strongly depends on A (or rg). Most 
remarkably, the asymptote D(oo) is systematically larger 
than the free diffusivity. Dp, for A < j/l/4, like for eccen- 
tric compartmentalized channels 43] . Note also that the 
slowly increasing branches of the curves plotted in Fig. 
|4]^b) tend to overlap and grow slower than for coaxial 
compartmentalized channels, where D(F) diverges like 

f2 gniis]. 

Like for channels I, two different transport mechanisms 
became apparent as we space the circular obstacles out 
along the channel axis, i.e., for xl » Ul- The diffusivity 
curve D{F) develops a secondary peak at lower F, which 
grows prominent on further increasing xl [Fig. [5jb)]. 
At variance with the high-F peak, the position of the 
emerging peak is rather insensitive io xl- Correspond- 
ingly, the mobility curves shift toward higher values, their 
most apparent feature remaining its non-monotonic F 
dependence. By closer inspection, for large x/,, in corre- 
spondence with the secondary D{F) peak one notices the 
appearance of a shoulder, or step, on the low-i^ raising 
branch of /x(F). 



IV. DEPINNING MECHANISMS 

The most remarkable property of transport in chan- 
nels I and II emerges for rectangular compartments with 
> Ul- As anticipated in Sees. |III A| and |III B[ the 
two-peaked structure of D{F) results from the compcti- 
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FIG. 5; (Color online) Driven transport in a channel II with 
yL = 1, -Do = 0.03, ro = 0.4 and different period xj^: particle 
mobility /i (a) and scaled difTusivity D/ Dq (b) vs. the scaled 
force F/ Dq. Dotted lines represent ^oo, Eq. Q, in (a) and 



Doo, Eq. 110 1, in (b) 



tion between the two depinning mechanisms illustrated 
in Fig. [6] The Brownian particle overcomes the obstacle 
thanks to the combined action of the drive F and the 
noise £,{t). The drive, in particular, on one side pushes 
the particle against the obstacles, thus enhancing the 
low-i^ pinning by the array (raising D(F) branch), on 
the other side guides the particle around the obstacles, 
so as to coUimate its trajectories through a coaxial [chan- 
nels I, D{oo) = Dq] or eccentric sequences of bottlenecks 
[channels II, D{oo) > Dq]. 

For relatively low drives, the particle circumvents a sin- 
gle circular obstacle, irrespective of the array geometry 
and the drive orientation, by diffusing a transverse dis- 




FIG. 6: (Color online) Depinning mechanisms, (a) The Brow- 
nian particle circumvents a single obstacle by diffusing in the 
transverse direction (vertical arrows) and then drifting along 
the axis driven by F (horizontal arrow). (b),(c) Trajectory 
focusing along the channel lane results from the competing 
action of longitudinal drift (horizontal arrows) and transverse 
diffusion (vertical arrows) across the array. The transverse 
diffusion length is of the order of half the channel width in 
(b) and half the bottleneck width in (c). (d) Example of min- 
imum longitudinal mobility orientation in a rectangular 2D 
array, see Sec. |V] 



tance rp over the time it takes to drift the same distance 
ro [Fig. p[a)]. On equating the transverse diffusing time, 
T± = rQ/2Z?0i and the longitudinal drift time, ry = ro/F, 
we locate the low-F depinning threshold at 



F^'^/Do = 2/ 



(5) 



We recall that excess diffusion peaks are characteristic 
signatures of depinning thresholds ^ 05 . We stress that 
this estimate for F^^^ is expected to apply to both types 
of channels, I and II, and to any geometry and size of 
their compartments (i.e., to any 2D array). As shown 
below, in the case of square compartments. Figs. [2jb) 
andjljjb), low- and high-F peaks cannot be separated. 
The estimate of Eq. ^ closely locates the low-i^ peaks 
of all D{F) curves reported in Figs, [sjjb) andjsjb). Note 
that in the limit of narrow bottlenecks, A ^ 0, F'-^^/Da 
approaches 4/^^, irrespective of [Fig. [5jb)]. 

For large drives, channels I and II behave differently. 
The bottlenecks in channels I are coaxial [Fig. [6]jb)]. This 
means that the driven Brownian particle gets trapped 
within a couple of obstacles if, during the drift time it 
takes to cover the distance separating the center of a 
bottleneck from the nearest pair of obstacles, xl — tq, 
it diffuses a transverse distance of the order of half the 
width of the (reduced) channel, 2/l/4. Equating the cor- 
responding times, T± = y'j^/'i2DQ and r|| = {x^ — rQ)/F, 
we obtain an estimate for the position of the high-F dif- 
fusivity peaks in channels I, namely. 



(6) 
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In channels II with A < ?/l/4, i.e. ro > Vl/^, the 
focusing action of the drive forces the particle through 
a meandering path. This action becomes effective when 
the particle drifts the distance separating a bottleneck 
from the opposing obstacle [Fig. |6|^c)], 3:^/2 — ro, while 
diffusing across the lane passing through the bottleneck 
of a (reduced) channel, i.e., a transverse distance of the 
order of half its width, A/2. The corresponding depin- 
ning threshold is thus estimated to be 



F(2.")/i?o = 4(xi-2ro)/A2 



(7) 



Contrary to F'^^\ the thresholds F^^'^) and F^^'") de- 
pend on the channel geometry, being different for chan- 
nels I and II, and linearly shift to higher values with 
increasing the length of the channel compartments, x/,, 
or shrinking the obstacle size, rg. Despite the rough es- 
timate of the characteristic transverse diffusing length, 
the thresholds of Eqs. Q and Q locate quite closely 



the large drive diffusion peaks in panels (b) of Figs. 2|5^ 
For square channel compartments, xl — Ul, with narrow 
bottlenecks, tq ~ yL/'2, the two depinning mechanisms 
are not clearly distinguishable, being the size of the ob- 
stacle and the array lattice constant of the same order. 
Correspondingly, the diffusivity curves exhibit only one 
peak. 

The limiting values of n{F) and D{F) for _F = and 
F — > oo are also of some interest. For F — 0, mobil- 
ity and diffusivity are related through Einstein's iden- 
tity, Eq. |4] Analytical expression for fj,Q in restricted ge- 
ometries can be obtained only under special conditions 
[551 For instance, in channels I and II correspond- 
ing to square arrays with narrow bottlenecks, A — > 0, /xq 
can be expressed in terms of the mean first exit time for 
the Brownian particle to escape the interstitial region de- 
limited by four nearest neighboring obstacles [Fig. ^h), 
inset]. Our numerical simulations show that, as expected, 
the mean exit time through all four openings, T4, is half 
the mean exit time through any pair of opposite openings, 
T2, with both T2 and decaying inversely proportional 
to the pore width, i.e., like 1/y/A. These results are well 
fitted by the law. 



T2 = 2t4 



(8) 



which is consistent with a recent analytical prediction 
[see Eq. (13) of Ref. @S]]. Moreover, the mean sojourn 
time in a square channel I compartment is 2t2 and, ac- 
cordingly, /io = xj^/4:DQT2, also in close agreement with 
the numerical data for small A plotted in Fig. [2j This ar- 
gument can be easily generalized to rectangular arrays, 
as long as the distance between neighboring obstacles 
remains sufficiently small compared to the lattice con- 
stants. 

The asymptotes ^(00) and D{oo) in channels II with 
^ < 2/l/4 also deserve attention. In the limit F — >■ 00, 
the pore eccentricity, which is responsible for the negative 
differential mobility of channels II, is eventually super- 
seded by the funneling effect due to the rounded shape 



of the obstacles. In such limit, the particle trajectories 
consist of straight parallel segments, which a free par- 
ticle traverses with speed F, and arcs of the obstacle 
boundaries, as shown in Fig. [ijc), where the driven par- 
ticle slows down. On averaging out noise fluctuations, 
the total time taken by the particle to cross one channel 
compartment, r^, can be derived analytically, so that, 
Moo = xl/tcF, that is 



XL 



V2e 



^Ini 

XL e 



-72^^ 



(9) 



with e — 2ro/?/L and 1/2 < e < 1. The predicted values 
for ^00 are indicated by dotted lines in Figs, [ijja) and 
[sja). Note that in channels II for A > 1/^/4 and in all 
channels I the traversal time is Tc = Xl/F and, therefore, 
A^oo = 1- 

A full analytic calculation of D{oo) proved to be a 
cumbersome task. As in channels II each bottleneck is 
faced by a blocking obstacle, one would expect that for 
A < ro the diffusivity D{oo) scales like ro/A, namely 
the ratio of the obstacle diameter to the vertical distance 
between two neighboring obstacles. Furthermore, on re- 
garding the channel as a ID discretized structure with 
spatial constant a;^, one can formally express D{oo) as 
a;^/2f, where f is a certain diffusion time constant, which 
we know to be proportional to the compartment volume 
[45] . Hence, even without calculating f, one concludes 
that D(oo) is proportional to xl- Therefore, we tried to 
reproduce our numerical data for D[oo) by means of the 
heuristic law 



D{oo) 



XL 

2yL 



ro 
A' 



(10) 



which turns out to work well both in Fig. Qb) and in 
Fig. [Sj^b). Note that, in the case of larger xl, Fig. [sf^b), 
where Eq. (10 1 seems to work not as well, a conclu- 



sive estimate of the asymptote D{oo)/Dq would require 
simulating even higher F/Dq values, which in practice 
becomes very difficult to implement numerically because 
this regime would require exceedingly small time steps. 



V. PARALLEL AND LATERAL TRANSPORT 



We address now the practical situation jT8H23t l25l - f29] 
when the orientation of the drive can be varied at will 
with respect to the array principal axes. The angle 
in Fig. [ija) can thus be rotated between and 7r/2. 
For a square lattice, = Ly, this corresponds to con- 
tinuously switch from a channel I (</> = 0) to a channel 
II ((/) = 7r/4) and then back again to the initial chan- 
nel I geometry (0 = tt/2). For a rectangular 2D lattice, 
Lx 7^ Ly, the two channels I have different width, that 
is, Ly for (f> = and for = 7r/2, while no channel 

II is recovered. Indeed, as F is directed along a diagonal 
of the rectangular lattice cell, the corresponding reduced 
ID channel looks similar to channel II in Fig. [ijc), but 
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FIG. 7: (Color online) Angular effects, (a) Migration an- 
gle; (b) longitudinal, fi±, and lateral mobility, fi\f, and (c) 
scaled longitudinal, D_i_/Do, and lateral diffusivity, D^^/Dq, 
of a Brownian particle driven at an angle <j) in rectangular 
arrays with different Lx (see legends). All quantities are plot- 
ted vs. for Ly = 1, ro = 0.4, and F/Dq = 10^. In these 
simulations we set Do = 1. 



for a relative shift of the lower and upper obstacle rows. 
The transport properties of such diagonal channels are 
qualitatively the same as for channels II of Sec. |III B[ 

In Fig. [7] we display some significant simulation data 
for driven Brownian transport in square and rectangular 
arrays. To this purpose we generalized the definitions of 



Eqs. ^ and ^ to extract the mobility and the diffu- 
sivity in the direction parallel and orthogonal to F. The 
corresponding four transport quantifiers are denoted by 
D|| for the longitudinal transport and by fi±, D± for 
the lateral or transverse transport. 

The salient properties of directed transport in 2D ar- 
rays are easily reconciled to the geometric properties of 
the particular system at hand. With focus on the (j) de- 
pendence of the transport quantifiers, one notices that: 

(i) For square arrays and vanishingly small drives, 
D±{0) = -011(0) and n±{0) — f^\\{0), irrespective of (f>. 
This is consistent with the D4 symmetry of the array 
lattice and with the fact that the mobility and diffusiv- 
ity components are defined in the long-time limit. Ap- 
plying a finite drive to a square array breaks the D4 
symmetry and leads to weaker angular symmetry re- 
lations, -D||,_l('0) = f^\\{-ip) = MlK-V"), and 
fi±{ip) — — "0), where ip = (j) — -k/A. Note that the 
lateral mobility is an odd function of tjj, all other trans- 
port quantifiers being even functions. 
(ii) Driven longitudinal transport in a square array is the 
least efficient for two special orientations, ±1/)*, or ((>* and 
7r/2 — (/)*, as signaled by the occurrence of two symmet- 
ric dips in and peaks in and -D_l [Figs. [7|b) and 
(c)]. The lateral mobility, vanishes for (f) = 0,7r/2 
(channel I geometry) and in correspondence with the dif- 
fusion peaks, -0 — zLtp* . The interpretation of such an- 
gular dependence follows from the argument leading to 
Eq. ([9]). A noiseless driven trajectory consists of straight 
segments parallel to F and circular arcs running around 
the obstacles. Each straight segment is tangent to the 
circumvented obstacle and impinges on the next obstacle 
that stands on the particle's way. The blocking action 
of this second obstacle is the largest, i.e., /iy has a min- 
imum, when the incident segment runs along one of its 
diameters [Fig. [6]jd)]. Under these circumstances, the 
trajectory can run around the blocking obstacle on ei- 
ther side with equal probability, i.e., fi± = 0, and the 
dispersion of the transport current in both directions, 
i.e. Dy and D±, is maximum. For narrow bottlenecks a 
simple calculation yields (f>* — arccos -^/l — {tq/ LyY, in 
agreement with the plots of Fig. [7] 

( Hi ) The lateral diffusivity in driven arrays is suppressed 
for increasing F (not shown). This effect signals that 
the particle is effectively channeled by oriented rows of 
bottlenecks [possibly at an angle with F, see Fig. [7|a)]. 
Such mechanism sets on when the driven particle crosses 
an array unit cell in a time much shorter than the time it 
takes to diffuse across a single bottleneck. The effective 
migration angle [35], a, of a particle driven across a 2D 
array does not necessarily coincide with the drive orien- 
tation. In fact, a = 0-|-arctan(/ij_//iy). This means that 
a and <f> coincide only for /ij^ — 0, that is for ?/; — 
and ±7r/4. For a square lattice the angular dependence 
of a(^) is plotted in Fig. J^a). Due to the angular sym- 
metry of ^J,± and /iy, it follows that a(V') + a(— V") = ^/2. 
The step-like structure of a(0) is controlled by the com- 
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mensuration of the lattice constants, as first reported in 
Ref. [28j. 

(iv) For rectangular arrays, the angular symmetries of the 
transport quantifiers in (i)-(iii) are broken. The nodes of 
the curves fi± {(j>) plotted in Fig. [T^b) increase in number 
with the ratio L^/Ly. Indeed, reflecting the lower rota- 
tional lattice symmetry, there exists more distinct critical 
angles, (/>*, that satisfy the maximal blocking condition 
of item (ii). However, as grows much longer than Ly, 
the lateral mobility for </) < 7r/4 gets suppressed, because 
channeling in the horizontal direction becomes more ef- 
fective than in the vertical direction. Correspondingly, 
the migration angle a{4>) develops more steps. For com- 
mensurate lattice constants, i.e., L^/Ly = n with n an 
integer, we counted exactly n steps of the a{(j)) curves 
within the interval < < 7r/4, each corresponding to 
a node of For it/A < (j> < tt/2, the step and node 

structure are smeared out when increasing the radius of 
the obstacles. Eventually, in the limit ^ Ly, lateral 
transport is appreciable only for (f> > tt/2 — (p*. Under 
such drive conditions, the curve boils down to a 

brad peak, which is seemingly independent of L^.. Cor- 
respondingly, the longitudinal mobility, /xy , is suppressed 
and the longitudinal diffusivity grows orders of magni- 
tude larger than Dq. 

VI. SUMMARY 

Forced transport across a 2D array was investigated for 
a simplified model where a single overdamped Brownian 
particle of negligible size is suspended in an unmovable 
interstitial fluid at fixed temperature. The particle is free 
to diffuse in the connected space delimited by circular 
reflecting obstacles of finite radius, arranged in a rect- 
angular lattice. The particle is only subject to thermal 
fluctuations and a homogeneous constant driving force. 

Brownian particle transport in the stationary regime 
was investigated by analyzing the dependence of the par- 
ticle mobility and diffusivity on the external drive (mag- 
nitude and orientation) and the array geometry (obstacle 
radius and lattice constants). Novel transport properties, 
including negative differential mobility, excess diffusion 



peaks, and unconventional asymptotic behaviors, which 
went unnoticed in earlier reports, have been detected by 
means of extensive numerical simulations. Such proper- 
ties have been explained in terms of two distinct depin- 
ning mechanisms: (i) trapping by a single obstacle, and 
(ii) correlated collisions against obstacle rows at an angle 
with the external force. The corresponding length-scales 
have been identified to be, respectively, the obstacle ra- 
dius (independent of the array geometry) and the effec- 
tive obstacle spacing (resulting from the combination of 
array geometry and drive orientation). 

All other effects, including asymmetry of the obstacles 
and of their spatial arrangement, particle-particle and 
particle-obstacle interactions, hydrodynamic corrections, 
chaotic and inertial dynamical terms, particle size and 
shape, have been ignored in this report. The generaliza- 
tion of our analysis to incorporate inertia and hydrody- 
namic corrections is the next step of this research project. 

Despite the various approximations and simplifications 
discussed above, the effects investigated here are robust 
enough to challenge experimenters investigating the dif- 
fusion of extended objects [lDl[5Tl|lSl|iai3Il[31l^or 
even point-like charge carriers [JTj in 2D arrays. More- 
over, the interplay of different diffusion length scales has 
been recently invoked also to explain certain features of 
the long-time self-diffusion of spherical tracer particles in 
periodic porous nanostructures [15] . 
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